typeset using JPSJ.sty <ver.l.Ob> 



00 
G\ 
OS 

Oh- 
(N 
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The temperature dependence of the correlation length, susceptibilities and the magnetic struc- 
ture factor of the two-dimensional spin-1 square lattice quantum Heisenberg antiferromagnet 
are computed by the quantum Monte Carlo loop algorithm (QMC). In the experimentally rel- 
evant temperature regime the theoretically predicted asymptotic low temperature behavior is 
found to be not valid. The QMC results however, agree reasonably well with the experimental 
measurements of La2Ni04 even without considering anisotropies in the exchange interactions. 
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Recently, field theoretical predictionsoErB 1 concern- 
ing the correlation length of the square lattice quan- 
tum Heisenberg antiferromagnet (QHA-J-were directly 
checked by experimental measurementsu'B^ .-and several 
quantum Monte Carlo (QMC) simulations BO &EP While 
in the case of spin S = 1/2 the validity of the predic- 
tions seemed supported by experiments, in the case £>f 
S = 1, experimental measurements for both La 2 NiOjj' 
and K^NiFjP turned out to be inconsistent with the the- 
oretical predictions. 

The inconsistencies were explained by noting that the 
theoretical low temperature expression is valid for the 
temperatures of the S — 1/2 experiments, but not valid 
in the temperature regime of the S = 1 experiments. 
While a number of simulations have been performed 
for S — 1/2. only a high temperature series expansion 
calculation^ .and. an effective high temperature theory 
(PQSCHAJli BEJ are available for S = 1 in the exper- 
imentally relevant temperature range. In this paper, we 
compute the correlation length and other thermal aver- 
ages for the S = 1 QHA on a square—Lattice using the 
quantum Monte. Carlo loop algorithmtir generalized to 
larger spins. E3E3'EIP This algorithm was implemented in 
continuous imaginary-time representation^ to eliminate 
the systematic error due to Suzuki- Trotter discretization 
of path integrals. 

When implemented in continuous imaginary time, the 
probability for the graph assignment in the cluster algo- 
rithms for larger spins becomes_miich simpler than the 
original discrete time version Jl3'E£P although the idea 
is essentially the same. As in the case of discrete time, 
we first extend the Hilbert space by expressing each spin 
operator by a sum of 25 Pauli spins: 
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where N is the total number of original spins. Since there 
are unphysical states in the new Hilbert space in which 
some of the spins have magnitude less than S, we must 
eliminate such states by applying projection operators. 
Here we use a representation where z-spin components 
are diagonalized. Our procedure for the graph assign- 
ment is as follows. For each pair of neighboring world 
lines and for each uninterrupted time interval during 
which spins on these world lines are antiparallel, we gen- 
erate "cuts" of worldlines with probability density J/2. 
At each cut, we reconnect pairwise the four end points 
created by the cut by two horizontal segments (Fig. 1). 
The application of the above-mentioned projection op- 
erator for a site is realized by choosing an appropriate 
boundary condition in the temporal direction. If the spin 
values (ofj ■, (r)) are the same at the four end points, (il) 
and (i2) at r = (3 and at r = 0, we choose a straight con- 
nection (connecting (il) at r = /3 to (il) at r = 0) and 
a cross connection (connecting (il) at r = f3 to (i2) at 
r = 0) with equal probability. Otherwise, we choose the 
unique method for connecting these four points pairwise 
so that the spin value at each connection point is contin- 
uous. Thus, we form many loops which are to be flipped 
with probability 1/2. This algorithm turns out to be 
more efficient than its discrete version. 

Using chiral|-perturbation theory, Hasenfratz and Nie- 
dermayer(HN)B) obtained the temperature dependence 
of the correlation length up to two-loop order for an ar- 
bitrary magnitude of spin: 
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We therefore consider (2SN) vertical lines along the 
imaginary time axis, each specified by two indices (?/Lt), 



8 2np s 

where t = k^T/2nps. The dependence on the magnitude 
of the spin S is only implicit through the «S-dependence 
of the the spin-stiffness constant p s , and the spin waise 
velocity c. For S — 1/2, the spin wave theory (SWT)ll$ 
values for p s and c were observed to be close to the QMC 
estimates. As SWT is better for larger spins, we can 
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Fig. 1. A cut and two horizontal segments. The two types of 
boundary condition in temporal direction are also illustrated. 



confidently use the SWT values for S = 1 in the following 
with 2irp s — 5.461J and he = 3. 067 J. These values 
agree with the result of a series expansion around the 
Ising limit, within statistical errors smaller than 1%. 

Equation (0) is valid in the renormalized classical 
regime where 



t<l. 



(3) 



An additional constraint on the temperature comes from 
the cutoff of the quantum fluctuations in the effective 
field theory once the extension in the imaginary time di- 
rection (3 becomes smaller than the lattice spacing a/Tic: 
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The approximation in the last term on the right hand 
side is again the leading order SWT result. 

In the case of 5 = 1/2 QHA, it was found by QMC sim- 
ulations, that eq. (g) is valid only at very large correlation 
lengthg§EP of the order of 100 lattice spacings or larger. 
In the experimentally relevant temperature regime the 
deviations, while clearly visible in the QMC simulations, 
are however smaller than the experimental errors. Thus, 
the theory and experiment agree for 5 = 1/2. nn 

The large discrepancies observed for 5 — laW are 
somewhat counterintuitive, since the theory based on a 
spin- wave picture should be better for larger spins. How- 
ever, this is not very surprising since eq. (g) reflects the 
quantum nature of the system and therefore eq. (|2J) ap- 
plies only if eq. (|4|) is satisfied. The region determined 
by eq. (|J) should be smaller for systems closer to the 
classical system. Therefore, since larger spins are more 
classical as noted previously in refs. [j], g, |[[l(j, the va- 
lidity of eq. (||) is restricted to even larger correlation 
lengths than for 5 = 1/2, much larger than accessible in 
experiments. 

As the analytic low temperature form is not valid in 
the experimentally accessible temperature regime, we 
performed QMC simulations over the temperature range 



t > 0.l8JkBT > J) in order to compare the experimen- 
tal datao' with the 5 = 1 QHA. The second moment 
correlation length £j, on a finite system of size L was 
determined from the static structure factor Sl(q) in the 
vicinity of Q = (n, 7r).EP The £l is calculated as follows. 

C L 2 = f(2n/L, 0) + /(0, 2n/L) - f(2ir/L, 2ir/L), (5) 
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This estimator for the correlation length (of a finite sys- 
tem) should be correct up to the fourth order in 2n/L. 
We use improved estimators to reduce statical errors. For 
a given temperature T, when the lattice size L is suffi- 
ciently large, £l converges to a size-independent value 
£(T), which we regard as the infinite-size limit. We find 
that the convergence is achieved to the accuracy deter- 
mined by the present statistical error when the condition 
L > 7£l is satisfied. All the following results are ob- 
tained under this condition. For each simulation we have 
performed 10 6 sweeps, after 10 4 thermalization sweeps. 

A selection of our results for the correlation length £, 
the staggered structure factor 5(7r, w) and the uniform 
susceptibility x are summarized in Table I. In Fig. 2, we 
plot our QMC results for thfi correlation length together 
with the experimental dataeP and theoretical predictions 
based on eq. (|2j). Our present estimates are in rough 
agreement with experimental measurements. Greven et 
alH and Nakajima et alJa have additionally proposed 
to include the effects of a small Ising anisotropy using 
a mean-field type correction to the theoretical isotropic 
results. This correction (eq. (7) of ref. Q) however makes 
the agreement worse, as can be seen from Fig. 2. 

Compared to theoretical predictions we find that the 
data deviate strongly from the low temperature for- 
mula_nf-.e£L (^). The effective PQSCHA approxima- 
tion JiSL3' 113 which is an effective high temperature the- 



ory, however agrees w^ 
ter than for S = 1/2 
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he QMC results much bet- 



Another violation of the theoretical predictionsa'LP was 
observed in the peak value of the staggered stijucture fac- 
tor S(tt,tt), which according to the theorjo tr c3 should 
scale as 



S{n, tt) = A2irM 2 £, 2 t 2 (l + Ct). 



(7) 



Here M is the staggered magnetization of the ground 
state and A and C are universal constants. For spin-1/2 
the leading t 2 form was confirmed by high temperature 
series. El Recent QMC datafl' at lower temperatures fit 
the above form very well in the temperature range t < 3 
with A « 4.0(1) and C w 0.5(1). However, experiments 
for both S =112, and 5 = 1 were better described by an 



empirical lav 



5(71", 71") oc £ 2 



(8) 



over the same temperature range. 

We applied a \ 2 analysis to check the consistency of 
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Table I. Correlation length §, magnetic structure factor S(n,n) 
and uniform susceptibility x as a function of temperature t = 
k B T/277p 3 . 



J/k B T 



t 



t 



S(n,n) 



0.10 


1.83 


20 


0.298(5) 


0.8915(1) 


0.051579(6) 


0.15 


1.22 


20 


0.396(5) 


1.0489(2) 


0.068503(9) 


0.20 


0.92 


20 


0.497(5) 


1.2481(3) 


0.08118(1) 


0.25 


0.73 


20 


0.610(5) 


1.5029(5) 


0.09038(1) 


0.30 


0.61 


20 


0.739(5) 


1.8342(7) 


0.09687(2) 


0.35 


0.52 


20 


0.894(6) 


2.272(1) 


0.10108(2) 


0.40 


0.46 


20 


1.079(6) 


2.854(2) 


0.10343(2) 


0.45 


0.41 


20 


1.307(6) 


3.648(2) 


0.10435(3) 


0.50 


0.37 


30 


1.604(3) 


4.776(1) 


0.104000(8) 


0.55 


0.33 


30 


1.982(8) 


6.392(4) 


0.10275(2) 


0.60 


0.31 


10 


2.482(5) 


8.795(2) 


0.100756(8) 


0.65 


0.28 


10 


3.14(1) 


12.46(1) 


0.09829(3) 


0.70 


0.26 


50 


4.06(1) 


18.25(2) 


0.09560(2) 


0.75 


0.24 


50 


5.30(2) 


27.60(3) 


0.09275(3) 


0.80 


0.23 


60 


7.06(2) 


43.30(5) 


0.09000(3) 


0.85 


0.22 


80 


9.52(3) 


70.03(8) 


0.08746(3) 


0.90 


0.20 


120 


12.98(4) 


116.4(1) 


0.08516(3) 


0.95 


0.19 


110 


17.97(5) 


198.8(3) 


0.08309(3) 


1.00 


0.18 


200 


24.94(7) 


344.0(4) 


0.08133(2) 
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o La 2 Ni0 4 measurements 

□ QMC data + Ising correction 

HN prediction 

PQSCHA approximation 



4 1/t = 2;ip s /k B T 5 



Fig. 2. The correlation length J as a function of t = k-QT/2Tvp s . 
Filled squares represent the results of «ur simulation, the open 
circles are experimental measurements. q) While the experimen- 
tal measurements agree roughly with our QMC results, they are 
incompatible with the QMC results corrected for a small Ising 
anisotropy (open squares). Compared to analytic calculations, 
we find that in this temperature range the low-temperature pre- 
dictions of eq. | (bL t tbyhed line) are not valid. The PQSCHA 
approximation J 1 Tl 12 [ F^P however, agrees well with both our data 
and the experimental measurements. 



the spin-1 data. Good fits were obtained for t < 0.28, 
with x 2 ~ 1 when we allow C to vary, and \ 2 ~ 2 with a 
fixed C = 0.5. The universal constant A was determined 
to be A — 4.1(1) and A — 4.5(1) , similar to the values 
obtained for S = 1/2. The discrepancies between the 
fits are non-universal effects caused by the rather high 
temperatures. 



Comparison of our data with the experiments is shown 
in Fig. 3 and we can see that for low temperatures the 
experimental data are consistent with the QMC results 
and eq. (M). The discrepancies that lead refs. H and 
to predict eq. (0) occur at higher temperatures where 
the experimental data has large error bars. In view of 
the precision of our QMC results and the large errors of 
the experimental results we suspect that, contrary to the 
suggestion of refs. H and 0, the deviations from eq. (Rl) are 
due to uncertainties in the experimental measurements. 
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Fig. 3. The ratio of structure factor peak value and square of the 
correlation length S(n, 7r)/§ 2 as a function of temperature. Solid 
circles represent the QMC measurements and open squares the 
experimental measurements. They agree at low temperatures, 
but differ at higher temperatures. The solid line is a fit of the 
QMC data to the theoretical low-temperature prediction. 



Finally, we present in Fig. 4 the uniform susceptibility 
for S = 1/2 together with previously published S = 1/2 
results. First we note that, as expected, the asymptotic 
low temperature behavior of the renormalized classical 
regime 



X 



3 V he ) 



(1 + t + t 2 ) 



(9) 



sets in at lower temperatures for S — 1 than for S = 1/2. 
It was discovered that for S = 1/2, the uniform suscep- 
tibility is the only quantitv for which a clear crossover to 
quantum critical behaviortf can be observed at interme- 
diate temperatures t ~ 1/3. The uniform susceptibility 
in the quantum critical regime is 



X = XI- + 



bt( 9 -^\ 
v he ) 



(10) 



For S = 1 



with a universal slope B s=s 0.26(1).&B 
however, as discussed above, non-universal corrections 
become important at lower temperatures. No quantum 
critical behavior was thus expected for S = 1. We can 
see in Fig. 4 that the uniform susceptibility for S = 1 
deviates from its universal quantum critical behavior at 
intermediate temperatures. Its slope is however still sur- 
prisingly close to the quantum critical one, indicating 
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that the non-universal corrections are still not very large 
for 5 = 1. 



1.5 



5 1.0 






0.5 



• spin-1 


/ o 


spin-1/2 


/ ° 
/ 




/ /° " quantum critical 




/ / 




■ f 


/; 


• 


m / ■:/ / 


. 


//i 


• 


fjf 


renormalized classical 


Jg^/ 


■ 


*gw ' 




j#tf&^ / 




'stir / 


■ 


/ 




/ 




/ 





0.0 



0.5 



1.0 



t = k B T/27ip s 



Fig. 4. The uniform susceptibility X as a function of temperature 
t for both spin 5=1 and spin S = 1/2 (taken from ref. ph. The 
solid line is the predicted low temperature form of the renor- 
malized classical regime. The dashed lines have the universal 
slope expected for the quantum critical regime. As expected, 
contrary to S = 1/2, no extended quantum critical regime exists 
for 5 = 1. However the slope is close to the quantum critical 
value as non-universal corrections are still small for S = 1. 



In summary, we have simulated the spin 5=1 quan- 
tum Heisenberg antiferromagnet on a square lattice in 
the experimentally relevant temperature regime k^T ~ 
J . We find a better agreement between the Heisenberg 
model and experimental data than is expected from the 
low temperature theory. However, in view of the exist- 
ing small discrepancies it may be necessary to perform 
simulations on a model with small anisotropies in the 
exchange interactions, and to critically check the data 
analysis of the experiments. 
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